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摘要 : 使 用 紫金 山 天 文 台 姚 安 观测 站 80 cm 望远镜 在 2022 年 2 月 至 4 月 对 矮 行 星 
Haumea 进行 了 两 轮 观 测 。 对 所 有 的 观测 图 像 使 用 干涉 条 纹 别 除 技术 ， 并 利用 相对 测 光 导出 
了 Haumea 的 仪器 星 等 ， 最 后 对 两 轮 光 度 测量 结果 进行 归 一 化 处 理 以 便 联合 分 析 。 通 过 相位 
弥散 最 小 化 (Phase Dispersion Minimization, PDM) 方 法 以 及 Lomb-Scargle 周期 图 法 分 别 求 得 了 
Haumea 的 自转 周期 。 两 种 方法 求 得 的 自转 周期 结果 仅 相差 0.072 s， 这 表明 所 求 周期 具有 良 
y= 好 的 一 致 性 。Haumea 光度 测量 结果 显示 明显 的 双 峰 自转 光 变 曲线 ， 自 转 周期 为 3.915 4 上 
= 0.0002 h， 峰 峰值 幅度 为 0.26+0.01 mag。 通 过 导出 的 自转 光 变 曲线 ， 最 终 估 计 了 Haumea 
由 于 自转 产生 的 光度 变化 对 位 置 测量 的 影响 最 大 为 -9~9 mas。 这 一 测 光 观测 对 后 续 高 精度 
天 体位 置 测 量具 有 基础 意义 。 
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太阳 系 柯 伊 伯 带 天 体 轨道 通常 在 海王 星 轨 道 以 外 ， 由 于 日 心 距 离 远 ,温度 低 ， 是 行星 吸 积 盘 原始 
遗迹 之 一 ， 因 此 携带 了 有 关 行 星 形成 的 物理 和 化 学 的 宝贵 信息 。 目 前 已 知 的 柯 伊 伯 带 天 体 超 过 4 000 
il (http ://www.minorplanetcenter.net) 。2004~2005 年 ，Brown HIBA (https ;//archive.org/details/howikille 
dplutow0Omike ) 和 Ortiz 41 BA ( http :// www. infoastro. com/200809/26ataecina-haumea. html ) 声称 发 现 了 矮 
> 行星 Haumea。 它 被 临时 命名 为 2003 EL61, J-F 2006 年 9 月 7 日 由 官方 小 行星 收录 ， 编 号 为 136108 。 
© Haumea ERAKHARRI —, HLÆR EE A Makemake (136472) 后 最 亮 的 海外 天 体 ， 其 
— 轨道 周期 为 284 地 球 年 ， 近 日 点 35 AU, WEAR 51 AU， 绝 对 星 等 0. 358 mag’! ， 半 长 轴 43.3 AU, 偏心 
率 0. 19， 轨 道 倾 角 28° 。 

光 变 曲线 是 获取 天 体 物 理性 质 的 主要 来 源 ， 为 小 天 体 及 行星 起 源 和 演化 的 研究 提供 了 重要 线索 。 
例如 , 文 [2] 推 导出 Haumea 自转 光 变 曲线 的 周期 为 3.915 4+ 0.000 2h， 峰 峰值 幅度 为 0. 28 + 0. 02 
mag， 颜 色 B-V=0.97+0.03mag。 文 [3] 使 用 夏威夷 大 学 2. 2m 望远镜 观测 ， 观 测 时 使 用 Kron-Cousins 
测 光 系 统 的 B，V，R 和 工 滤 光 片 ， 导 出 周期 为 3. 915 5+ 0. 000 1 h 的 双 峰 自转 光 变 曲线 ， 其 峰 峰 值 幅 
EX 0.29 + 0. 02 mag, 自转 周期 短 ，Haumea E JEFE HG ERA (Jacobi ellipsoid) ， 假 设 其 已 经 达到 流体 
静 力 平衡 ， 文 [2-3] 结 合 这 两 个 特性 推测 其 质量 密度 p~2 500 kg + m”, 

2005 年 ， 文 [14] 发 现 了 Haumea 的 两 颗 卫 星 Hi'iaka 和 Namaka， 这 两 颗 卫 星 的 星 等 较 Haumea 分 
别 暗 2. 98 +0. 03 mag 和 4.6mag。 由 于 地 面 望远镜 观测 受到 大 气 的 影响 ， 观 测 时 通常 不 能 分 状 Haumea 
及 其 卫星 ,但 是 天 体 测量 结果 可 以 用 于 确定 光 心 的 轨道 。 根 据 文 [5] ， 光 心 轨道 是 双星 系统 相对 轨道 
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的 缩放 版 本 ， 因 此 可 以 通过 光 心 的 运动 轨道 初步 推断 双星 系统 的 相对 轨道 。 于 是 ， 文 [5] 结 合 径 向 速 
度 测 量 结果 、 历 史 测量 结果 以 及 其 他 所 有 可 用 的 天 体 测量 结果 (包括 相对 位 置 测 量 值 、Hipparcos 中 间 
数据 以 及 追踪 光 心 移动 的 测 光 观测 结果 ) ， 对 双星 系统 GJ67AB 进行 轨道 和 质量 求解 ， 但 是 质量 测量 
的 精度 仍然 主要 受到 相对 位 置 测 量 的 影响 。 然 而 ， 对 于 光 变 幅度 较 大 的 天 体 ， 位 置 测 量 受 到 光度 变化 
的 影响 。 参 考 文 [5] 的 工作 ， 在 求解 Haumea 与 卫星 的 光 心 轨道 时 也 会 受到 相同 的 影响 ， 我 们 需要 先 
估计 影响 范围 。 

基于 紫金 山 天 文 台 姚 安 观 测 站 80 cm 望远镜 观测 的 Haumea CCD 图 像 资 料 ， 本 文 推导 Haumea 的 
光度 变化 ， 为 后 续 改进 相对 位 置 测 量 英 定 基础 。 


1 观测 资料 和 数据 处 理 


我 们 使 用 紫金 山 天 文 台 姚 安 观测 站 80 cm 望远镜 (详细 参数 如 表 1) AREF BE Haumea 进行 两 轮 测 
光 观 测 A 观测 时 使 用 Johnson-Cousins I EEF , 具体 的 观测 资料 情况 如 表 2。 
我 们 利用 Maxim DL ( https ://diffractionlimited. 


com/maxim-dl/) 软 件 对 获得 的 原始 图 像 进 行 预 处 表 1 望远镜 及 CCD 详细 参数 

理 主要 包括 本 底 i Re Hy AE Hye TEL ` 由 于 反 Table 1 Specifications of the telescope and its CCD 
射 光 在 CCD 硅 片 薄 层 间 来 回 反射 ， 形 成 干涉 条 纹 ， Items Parameters 
导致 图 像 背 景 不 均匀 。 尤 其 对 于 Johnson-Cousins I Focal length 800 cm 

滤 光 片 ，CCD 硅 片 的 吸收 率 随 着 波长 增加 迅速 下 人 pices 

降 ， 有 相当 一 部 分 低能 光子 没有 被 完全 吸收 之 前 ， gy i 

在 CCD 硅 片 的 薄 层 来 回 反 射 好 几 次 ， 产 生 干 涉 条 

纹 ， 不 利于 精确 的 光度 测量 。 因 此 需要 对 拍摄 的 图 eee 
像 进行 处 理 ， 剔 除 干 涉 条 纹 ， 具 体 的 步骤 如 下 。 首 Size of pixel 13. 5 wm x 13. 5 pm 
先 ， 从 每 晚 观测 的 CCD 图 像 中 选择 1 或 2 幅 图 像 。 Size of CCD array 2 048 x 2 048 
接着 ， 对 筛选 的 CCD 图 像 (剔除 了 30 ) 使 用 “3-2” pr ee 0. 348"/pixel 


公式 (3 倍 中 值 -2 倍 均值 ) 进行 天 空 背景 估计 ， 得 
到 每 幅 图 像 的 背景 估计 值 。 所 有 筛选 的 CCD 图 像 EEEE 

减 去 各 自 的 背景 估计 值 后 ， 将 这 些 图 像 每 一 个 对 应 E pannel ee 
REACHES, BGM Ee eee 
纹 模板 对 应 像素 的 灰 度 值 。 对 于 构造 的 干涉 条 纹 模 a p ne 
板 ， 需 要 乘 上 一 个 比例 因子 以 获得 归 一 化 干涉 条 纹 wu) ous me 
模板 ， 该 比例 因子 为 360 s 除 以 待 去 除 干涉 条 纹 图 2022-02-28~ 31 Johnson-Cousins I 360 
像 的 曝光 时 间 。 最 后 ， 从 待 去 除 干涉 条 纹 的 原始 图 70770-04 

像 中 减 去 归 一 化 干涉 条 纹 模板 ， 得 到 剔除 干涉 条 纹 2022-04-21~ 82 Johnson-Cousins 1 360 
后 的 图 像 。 典 型 的 处 理 结果 如 图 1。 2o 


2 光度 测量 结果 与 分 析 


2.1 光度 测量 数据 

本 文 使 用 相对 测 光 技术 获取 目标 的 相对 仪器 星 等 。 由 于 两 轮 观 测 时 视 场 中 的 星 像 较 少 ， 而 且 在 每 
轮 观 测 期 间 ， 目 标 仅 运 动 了 一 小 段 距 离 ， 因 此 ， 我 们 可 以 使 用 相同 的 参考 星 进行 测 光 :7 。 具 体 地 ， 
在 第 1 轮 观测 时 ， 选 择 3 颗 参 考 星 ， 在 第 2 轮 观测 时 ， 选 择 5 颗 参 考 星 。 所 有 的 参考 星 都 是 恒星 且 光 
度 测量 稳定 。 两 轮 观测 中 典型 的 目标 和 参考 星 的 图 像 如 图 2， 其 中 ， 目 标 Haumea 使 用 绿色 圆圈 标记 ， 
参考 星 使 用 红色 圆圈 标记 。 
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<5 Haumea 


(a) (b) (c) 


图 1 (a) 存 在 明显 干涉 条 纹 的 原始 图 像 ;， (b) 和 干涉 条 纹 图 像 ; (IRTIRA g, At Haumea 
位 于 图 像 中 心 附近 
Fig.1 (a) A typical frame before the fringe-removal; (b) fringing pattern; (c) a defringed frame, and the target 


Haumea is located near the center of CCD frame 


a gi | (b) = 


图 2 观测 拍摄 的 典型 视 场 图 像 。(a) 第 1 轮 观测 获得 的 图 像 ; (b) 2 轮 观 测 获得 的 图 像 。 两 张 图 像 剔除 了 干涉 条 纹 
Fig.2 Example of fields of view obtained during our observations. (a) Obtained in the first run; (b) obtained in the second 


run. Fringes are removed from both frames 


我 们 使 用 二 维 高 斯 定 心 算 法 获取 目标 星 和 参考 星 的 像素 坐标 ， 采 用 “3-2” 公 式 对 目标 星 和 参考 
星 周 围 环 形 区 域 天 空 背景 进行 估计 ， 使 用 photutils( https :A/pypi.org/project/photutils/) 测 光 库 提供 的 孔 
径 测 光 程 序 对 目标 星 和 参考 星 进 行 光 度 测量 。photutils 孔径 测 光 程 序 允 许 我 们 使 用 一 系列 不 同 的 孔 
径 ， 选 择 光度 测量 精度 最 高 的 孔径 作为 最 优 孔 径 ， 得 到 每 个 夜晚 目标 星相 对 仪器 星 等 。 误 差 主 要 来 源 
于 孔径 的 目标 信 噪 比 ， 同 时 我 们 也 考虑 了 目标 与 参考 星之 间 因 颜色 不 同 带 来 的 误差 。 

由 于 没有 进行 标准 测 光 ， 为 了 将 两 轮 观 测 的 测 光 结果 联合 起 来 进行 周期 分 析 ， 根 据 文 [8] ， 我 们 
分 别 拟 合 两 轮 观 测 求 出 的 目标 相对 仪器 星 等 ， 拟 合 结果 之 间 存 在 一 个 常数 偏差 ， 最 后 使 用 这 个 常数 改 
正 求 出 的 目标 相对 仪器 星 等 。 每 次 拟 合 采用 一 轮 多 晚 观 测 的 数据 ， 而 不 是 每 个 夜晚 的 观测 数据 。 
2.2 自转 周期 

为 了 测定 准确 的 Haumea 自转 周期 ， 本 文 使 用 两 种 周期 求解 方法 ， 分 别 是 相位 弥散 最 小 化 方法 以 
及 Lomb-Scargle 周期 图 法 。 两 种 方法 求 出 目标 Haumea 自转 周期 结果 和 讨论 分 析 如 下 。 

相位 弥散 最 小 化 方法 由 文 19] 在 1978 年 提出 。 相 位 弥散 最 小 化 方法 通过 求 统 计量 9 的 最 小 值 获取 
周期 ,9 定义 为 


0= =， (1) 


3 期 李 
N M 

pCR 3 Dyas 

9 i=l 2_ j=l 


ao = N-1 a4 M 


, (2) 


(ORP o 为 观测 的 w 个 离散 数据 点 的 方差 。 
(2) 式 中 >x; 是 指 在 i 时刻 观测 的 目标 星 等 ,x 为 
N 个 数据 点 的 均值 ; 从 V 个 数据 点 中 选择 1 个 
有 着 相同 相位 D, 的 样本 ， 每 个 样本 中 包含 数据 
SIDS n; 8? 为 第 7 个 样本 的 方差 。 通 常 ， 
9 值 在 0 到 1 之 间 。 如 果 求 出 的 是 正确 的 周期 ， 
那么 0 会 达到 一 个 局 部 最 小 值 ， 接 近 0。 使 用 相 
位 弥散 最 小 化 方法 计算 得 到 的 周期 如 图 3。 在 频 
率 为 12. 259 39(cycles/day ) 处 达到 局 部 最 小 值 ， 
最 小 值 约 为 0. 17。 根 据 自 转 频 率 O 与 自转 周期 
尸 之 间 的 关系 Q2=24/P， 可 以 求 出 自转 周期 尸 = 
1.957 68h, 

Lomb-Scargle 周期 图 法 5 适用 于 测定 非 均 
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0 5 10 15 20 25 
Frequency / (cycles /day) 
图 3 相位 弥散 最 小 化 方法 计算 的 9， 即 相位 弥散 值 作 为 
频率 的 函数 。 在 频率 为 12. 2$9 39( cycles/day) 时 出 
现 了 一 个 局 部 最 小 值 
Fig.3 Dispersion, 0, calculated by the phase dispersion 
method as a function of frequency. The data shows one 


clear local minimum located at 12.259 39 cycles/day 


匀 采 样 数据 的 周期 。 该 方法 基于 传 里 叶 谱 分 析 ， 考 虑 非 均匀 采样 数据 ， 并 赋予 每 个 数据 点 权重 而 不 考 


虑 时 间 间 隔 。 假 设 有 N 个 数据 点 万 =h(4)(j = 0， 


1 [LÈ h -hcoso -nF [E C- A) sino(y, -7)] 
+ 


…, NN - 1) ,4 为 实际 测量 时 刻 ， 按 昭 


Poo) 2 


定义 Lomb 归 一 化 周期 图 ， 其 中 ，o 为 角 频 率 ， 
o 为 数据 的 方差 ,7 定义 为 
> sin2ot, 
È, ,cos2or, 
Lomb 归 一 化 功率 谱 最 大 值 对 应 的 值 为 求解 的 最 
佳 周 期 。 使 用 Lomb-Scargle 方法 计算 得 到 的 周期 
如 图 4。 图 4 展示 了 一 个 明显 的 最 大 峰值 在 自转 
频率 为 12. 259 40 + 0. 000 04( cycles/day) ， 可 以 
求 出 自转 周期 P=1.9577+0.0001h。 与 相位 弥 
散 最 小 化 方法 求 出 的 结果 仅 相 差 0.072 s， 表 明 
两 种 方法 求 出 的 结果 具有 良好 的 一 致 性 。 
2.3 自转 光 变 曲线 分 析 

根据 文 [12] ， 对 于 像 Haumea 这 种 因 快 速 旋 
转 而 变形 拉 长 成 椭 球 体 的 TNO ( Trans-Neptunian 
Object) ， 其 完整 的 光 变 曲线 应 包含 两 个 极 大 值 
和 两 个 极 小 值 ， 极 大 值 代表 天 体 的 前 后 两 个 侧 视 


tan(2@7T) = 


(4) 


T È , cos’; -T) 


(3) 


> ' sin w(t -T) 
J p 


100 - 


Lomb-Scargle power 


0 5 10 15 20 25 
Frequency / (cycles /day) 
图 4 Lomb-Scargle 功率 谱 。 图 中 显示 最 大 峰值 在 频率 为 
12. 259 40( cycles/day) 处 
Fig.4 | Lomb-Scargle power spectrum computed from the 
combined data. The data shows one clear peak located 


at 12. 259 40 cycles/day 


面 ， 而 极 小 值 是 两 端的 视图 。 因此 ， 将 2. 2 节 求 解 的 Haumea 自转 周期 拉 长 至 两 倍 ， 求 得 Haumea 的 
自转 周期 7=3. 915 4+ 0. 000 2 php， 该 结果 与 文 [2-3] 的 结果 吻合 。 
我 们 使 用 一 般 的 线性 最 小 二 乘法 拟 合 测 光 数据 (剔除 了 30) ， 拟 合 模型 为 
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RI 拟 合 模 型 系数 详细 说 明 


3 
27 27 
.二 C+ ‘sin(1— t.) +b. I — t. 
a 2 区 T ti) bcos (j T 1) | . ta) Table 3 Specifications of the coefficients 


FEA, 为 观测 时 刻 的 JD 时 间 ( 校 正 了 光 行 时 ) 。 各 系数 拟 合 of the fitting model 
结果 和 误差 如 表 3。 观 测 数据 点 和 拟 合 模型 光 变 曲线 如 图 5, Coefficients Value Error/mag 
图 中 显示 两 个 极 大 值 分 别 为 0. 12 mag 和 0. 10 mag， 两 个 极 小 c -2. 801 0. 002 
值 分 别 为 0. 14 mag 和 0. 09 mag ， 对 应 的 如 =252， 拟 合 结果 服 a, -0. 023 0. 002 
KA 由 度 为 205 的 x? 分 布 。 x? 定义 为 b, —0. 009 0. 003 
N Yı 一 m; 2 
Ys > | i | , (6) a> 0. 099 0. 003 
iL a; b, -0. 046 0. 002 
其 中 ,VN 为 拟 合 数据 点 总 个 数 ，y, 为 第 i 个 数据 点 对 应 的 相 A _0 006 0 002 
HYI EE. m> + yee, > ;个 数 握 占 的 测 i 
ETER m; 为 (5) 式 计算 的 结果 ; O; 为 第 1] 数据 点 的 测 b, —0. 002 0. 002 
量 误差 。 
= 2022-02-28~2022-03-04 Yaoan-80 cm 
F -0.2 F - 2022-04-21~2022-04-29 Yaoan-80 cm 
g f 
& oot 
£ 
2 0.1 上 
2 + 
E tT f 
2 02 上 
0.0 0.2 0.4 0.6 0.8 1.0 


Rotational phase 


图 5 Haumea 观测 光 变 曲线 和 拟 合 模型 光 变 曲线 。 散 点 为 观测 数据 ， 红 线 为 拟 合 模型 光 变 曲线 
Fig.5 The observational and fitting model lightcurves of Haumea. The scattered dots are observed data, 


and the red solid line are fitting lightcurve 


通过 图 5 可 以 看 出 ， 模 型 光 变 曲线 拟 合 了 大 部 分 的 观测 数据 点 。 少 部 分 数据 点 残 差 较 大 ， 可 能 是 
因为 时 间 间 隔 内 采样 稀 下 或 者 数据 点 本 身 误 差 较 大 ， 也 可 能 由 于 受到 Haumea 旁边 两 颗 卫 星 Hi'iaka 
和 Namaka 光 变 的 影响 ， 该 数据 点 弥散 较 大 。 


3 ” 光 变 对 位 置 测量 的 影响 


对 于 像 Haumea 这 种 光度 变化 幅度 较 大 的 天 体 ， 光 度 会 
对 位 置 测量 造成 一 定 的 影响 。 由 于 Namaka 较 暗 ， 参 考 文 
[1] 的 工作 ， 只 考虑 Haumea 与 卫星 Hi'iaka 构成 的 双星 系统 
(如 图 6) 。 在 天 空 平 面 上 ， 设 主星 与 卫星 之 间 的 角 距 离 为 
d; 质心 与 光 心 之 间 的 角 距 离 为 A; P AS 分 别 代表 主星 和 


卫星 ; Cy 和 Ci 分 别 是 质心 和 交心 。 根据 杠杆 原理 ， 
1 1 图 6 Haumea RZ AR‘! 
A= | | l 


13] 


LFA 一 Lr, (7) Fig.6 The binary system of Haumea! 


FUP, r =L,/L, sry =m,/m,; 二 和 到 分 别 是 主星 和 卫星 的 亮度 ; m, Mm, 分 别 是 主星 和 卫星 的 质量 。 
假设 卫星 Hi'iaka 的 光 变 幅度 不 变 ，A 的 变化 主要 由 Haumea 光 变 引起 ， 可 以 根据 误差 估计 公式 估计 4 


ae 5 ae 1 1 
受 光 变 影响 的 变化 范围 。 令 (8) 


= 9 
l+r 1l1+ry 


3 期 李 有 丹 等 面向 天 体 测量 的 矮 行 星 Haumea 相对 测 光 试验 

则 有 A = kd， 那么 A 受 光 度 的 变化 影响 为 A’=(k+6k)d, (9) 

ek re À L, ôL 

k 的 变化 只 与 主星 和 卫星 的 亮度 变化 有 关 ， 因 此 有 ks- a (10) 
> > 


其 中 ， Ly 是 主星 与 卫星 的 亮度 总 和 ; 根据 文 [4] ， 卫 星 Hi'iaka 的 亮度 是 主星 Hamuea 的 5.9%, Ath 
r,=16.9, L/Ly=1/17.9, 根据 文 [14] ，ry =223. 80。 根 据 星 等 转换 关系 式 有 6Ls/Ls=10” 4-1, 其 


H, Am 为 Haumea 光 变 曲 线 上 对 应 的 星 等 值 。 


从 IMCCE 历 表 ( https://www.imece.fr/) 可 以 7O00 U A 

获得 任意 时 刻 主 星 与 卫星 之 间 的 距离 4。 在 整个 -0.007 5 F ik /\ 
观测 期 间 ， 选 取 Haumea 与 卫星 Hi'iaka 的 距离 。 0.0050 ， i 4 j 
BE (d= 1.446") 的 时 间 节 点 ， 每 隔 1 min 取 数 ~-0.0025 F $ i ， i 
据点 ， 一 共 300 个 数据 点 , HERE T Haumea 一 个  Ș 00000| ` j 
自转 周期 (T= 3.915 4h), RIE 2.3 节 求 得 的 。 Š ooosk $ ; ， j 
Haumea 光 变 曲线 图 (图 5)， 我 们 可 以 获得 Haumea 0.005 0 上 VU i } 
光度 的 变化 对 位 置 测量 的 影响 曲线 图 (如 图 7) 。 wiet VV 

a 其 中 ， 横 坐标 为 Haumea 的 自转 相位 ， 与 图 5 的 gogo es i i | | 
横 坐 标 相 对 应 ， 纵 坐标 为 该 时 刻 因 光度 变化 ， 0.0 02 04 0.6 0.8 1.0 
Haumea 位 置 测 量 受 到 的 影响 。 可 以 看 出 ，Haumea Rotational phase 

= 光度 变化 对 位 置 测 量 的 影响 最 大 可 达到 -9~9 图 7 Haumea 光度 变化 对 位 置 测 量 的 影响 

= mas。 因 此 ， 利 用 Haumea 的 光 变 特征 改正 其 位 Fiet Tueee of Hannei lehtas aotik 


on position measurement 


置 测量 是 有 必要 的 。 


4 总 结 与 展望 


我 们 对 紫金 山 天 文 台 姚 安 观测 站 80 cm 望远镜 拍摄 的 矮 行 星 Haumea CCD 图 像 剔 除了 干涉 条 纹 ， 


> 使 用 相对 测 光 技术 获得 Haumea 的 相对 仪器 星 等 ， 通 过 联合 两 轮 观测 资料 探究 了 其 光 变 情况 。 本 文 使 
Se 用 相位 弥散 最 小 化 方法 以 及 Lomb-Scargle 周期 图 法 计算 Haumea 的 自转 周期 ， 结 果 显 示 一 条 明显 的 双 
© 峰 自转 光 变 曲线 ， 周 期 为 3.915 4+ 0. 0002h， 峰 峰值 幅度 为 0. 26+0. 01 mag; 最 终 估 计 了 Haumea 由 
= 于 自转 产生 的 光度 变化 对 位 置 测量 的 影响 最 大 为 -9~9 mas。 在 后 续 的 工作 中 ， 我 们 将 结合 Haumea 光 
= 度 测量 与 位 置 测量 ， 利 用 其 光度 变化 特征 改进 位 置 测 量 精 度 。 
致谢 : 感谢 紫金 山 天 文 台 姚 安 观测 站 80 cm 望远镜 运行 组 全 体 成 员 。 感 谢 郑 中 杰 在 本 工作 期 间 提 


供 的 帮助 。 
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Astrometry-Oriented Relative Photometry Test of Dwarf Planet Haumea 
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Abstract: The 80 cm telescope at Yaoan Station, Purple Mountain Observatory was used to observe the 


dwarf planet Haumea for two runs from February to April 2022. Interference fringes were removed from all the 


frames, and the instrument magnitude of Haumea was derived by relative photometry technology. Then, the 


photometric data were normalized for combining two observation runs to analysis. The rotation period of Haumea 


was obtained by the PDM (Phase Dispersion Minimization) method and the Lomb-Scargle Periodogram , 


respectively. The difference between the results by the two methods was 0.072 s, which showed good 


consistency. Our photometric results reveal a distinct, double-peaked rotational light curve with a period of 
3.915 4+ 0.000 2 h and peak-to-peak amplitude of 0. 26 + 0.01 mag. Through the derived rotational light 


curve, it is finally estimated that the maximum influence of the photometric change of Haumea due to rotation 


on the astrometry is -9~9 mas. The test has a basic significance for future high-precision astrometry. 


Key words: relative photometry; Kuiper belt; period determination; CCD image processing 


